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Abstract 

We review the status of — * u T flavor transitions of atmospheric neutrinos in the 92 kton-year 
data sample collected in the first phase of the Super-Kamiokande (SK) experiment, in combination 
with the recent spectral data from the KEK-to-Kamioka (K2K) accelerator experiment (including 
29 single-ring muon events). We consider a theoretical framework which embeds flavor oscillations 
plus hypothetical decoherence effects, and where both standard oscillations and pure decoherence 
represent limiting cases. It is found that standard oscillations provide the best description of 
the SK+K2K data, and that the associated mass-mixing parameters are determined at drier (and 
AT DF = i) as: Am 2 = (2.6 ± 0.4) x 10~ 3 eV 2 and sin 2 29 = LOOt^. As compared with standard 
oscillations, the case of pure decoherence is disfavored, although it cannot be ruled out yet. In the 
general case, additional decoherence effects in the v„ — ► v T channel do not improve the fit to the Sk 
and K2K data, and upper bounds can be placed on the associated decoherence parameter. Such 
indications, presently dominated by SK, could be strengthened by further K2K data, provided 
that the current spectral features are confirmed with higher statistics. A detailed description of 
the statistical analysis of SK and K2K data is also given, using the so-called "pull" approach to 
systematic uncertainties. 

PACS numbers: 14.60.Pq, 13.15.+g, 03.65.Yz, 95.55.Vj 
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I. INTRODUCTION 



In its first phase of operation (years 1996-2001), the Super-Kamiokande (SK) experiment 
has provided, among other important results, compelling evidence for atmospheric disap- 
pearance 0,0. This evidence, now firmly based on a high-statistics 92 kton-year exposure 
has not only been corroborated by consistent indications in the MACRO 5] and 
Soudan 2 [(J atmospheric neutrino experiments, but has also been independently checked 
by the first long-baseline KEK-to-Kamioka (K2K) accelerator experiment [7|, |8(, using SK 
as a target for z/^'s produced 250 km away with (E u ) ~ 1.3 GeV. 

Neutrino flavor oscillations, originated by nonzero mass-mixing parameters (Am 2 , sin 2 20) 
in the z/ M — > v T channel, provide by far the best and most natural (and probably unique) 
explanation for the observed disappearance [iNA- Among the "exotic" scenarios which 
have been proposed as radical alternatives (see^li , M EoT] for reviews), at present only 
hypothetical neutrino decoherence ED (see also BEl) appears as a possible "survivor" in 
the official SK data analysis [3. ll4lTl5l]. All the other models envisaged so far to challenge 
the standard picture are strongly rejected as dominant explanations, and are bound to give 
at most sub dominant effects (see, e.g., 0,0,0). 

In this paper, we review the phenomenological status of both standard oscillations and 
decoherence effects in the Vj, — > v T channel, in the light of the latest SK atmospheric zenith 
distributions of leptons |H 0, El and of the first spectral results from the K2K experiment 
[1, El- In Sec. II we briefly review the adopted theoretical framework and describe the data 
set used. In Sec. Ill and IV we show, respectively, the results of our analysis for the cases 
of standard oscillations and pure decoherence, while in Sec. V we consider the most general 
case (oscillation+decoherence). Conclusions are given in Sec. VI. 

We find that standard oscillations provide the best fit to the SK+K2K data, and that 
the mass-mixing parameters are determined at ±lcr as 

Am 2 /eV 2 = (2.6 ±0.4) x 10" 3 , (1) 
sin 2 26 = 1.001°;™ , (2) 

with errors scaling linearly, with good accuracy, up to ~3o~. Conversely, possible additional 
decoherence effects are generally disfavored, although not very strongly. We discuss how the 
upper bounds on such effects, currently dominated by SK, can potentially be improved by 
further K2K spectral data. The technical details of our SK and K2K data analysis are finally 
described in two appendices, using the so-called "pull" approach to systematic uncertainties. 



II. THEORETICAL FRAMEWORK AND EXPERIMENTAL DATA 

Oscillating neutrino systems may be thought as extremely sensitive interferometers with 
long arm lengths. When the interferometer beams interact with a background (e.g., matter), 
the interference pattern may be altered, and may even disappear if the background is "fuzzy." 
In theories of quantum gravity, it has long been speculated that the space-time itself may 
act as a "fuzzy" or "foamy" background for any propagating particle, leading to possible 
decoherence and dissipative phenomena on macroscopic scales (see, e.g., the bibliography in 
lH .) In this case, the neutrino evolution has to be described in the language of quantum 
open systems interacting with a generic environment 2, 13| . 

In Ref. |ll| the formalism of quantum open systems has been applied to the case of two- 
family Va — > v T transitions of atmospheric neutrinos. Assuming entropy increase and energy 



2 



conservation in the neutrino subsystem, the evolution equation was shown to depend on 
only one additional (decoherence) parameter, here denoted as /i 2 , in addition to the usual 
mass-mixing parameters (Am 2 , sin 2 29). 1 With the further assumption of Lorentz-invariance 



of n , it was shown that the Vn survival probability takes, in natural units, the form 



sin 2 29 



e 2E v cos 



'Am 2 !/ 
2E V 



(3) 



where E u (L) is the neutrino energy (pathlength) , and /i 2 has the dimensions of a squared 
mass. In the absence of first-principle calculations, fi 2 must be considered as a purely 
phenomenological parameter, to be constrained by experimental data. 

The above equation has two interesting limits. The first is reached for /i 2 
sponding to the well-known case of standard oscillations, 



0, corre- 



^7(Am 2 ,i 



sin 2 29 
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which perfectly fits atmospheric v data with Am 2 ~ 3 x 10 3 eV 2 and sin 2 29 < 
second limit is reached for /z 2 = 0, corresponding to the case of pure decoherence 



(4) 



1. The 



P, d , ec (/z 2 ,fl) 



1 - 



sin 2 29 







1 — exp ^ 


2E U J 



(5) 



where the oscillation pattern is completely absent and the disappearance rate becomes 
monotonic in L/E v . Surprisingly, the pure decoherence case turned out to fit well the 52 
kton-year data SK atmospheric data for /i 2 ~ 3 x 10~ 3 eV 2 and sin 2 29 ~ 1, with a x 2 
comparable to the standard oscillation case Only with almost doubled SK statistics 

(92 kton-year exposure y|_) some weak indications have started to emerge against the pure 
decoherence scenario [3, EH, Ell • 

It is important to note that neutrino decay ^1 can lead to an exponentially decreasing 



P w qualitatively similar to Eq. JSJ). The decay and decoherence scenarios cannot thus be 
easily distinguished in the disappearance channel. However, they can be distinguished 
through the appearance mode, i.e., through neutral current (NC) events and r appearance 
events in SK j3j. In fact, while in the decoherence case the total number of active neutrinos 
is conserved, in the decay scenario it decreases with LjE v . The nonobservation of NC event 
suppression and the indications for r appearance in SK can thus be used to reject decay 
|3|, y, ll5j but not decoherence effects. 

We remark that, in the general case [Eq . (13 1 ]. decoherence effects basically acts as a 
damping factor for the oscillating term Il8l | . Testing atmospheric and long-baseline 
accelerator neutrino data through Eq. (j3J) amounts thus to test how well the oscillatory 
pattern is favored by (or hidden in) the data. Therefore, the decoherence scenario represents 
a useful benchmark, although its theoretical motivations are admittedly weaker than those 
supporting standard oscillations. In practice, in fits with unconstrained (Am 2 , /i 2 , sin 2 29), 
the emergence of the standard oscillatory pattern should be signalled by a preference for 
small values of fi 2 . As we shall see, SK and K2K consistently favor fi 2 ~ 0; however, the 
absence of a clear oscillation pattern in the data makes this indication not very strong yet. 



The relation with the parameter /3 in 11] is \i l = 2/3 
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We conclude this section with a review of the SK and K2K data sets used for our analysis 
(see the appendices for more details). Concerning SK atmospheric neutrino data (92 kton- 
year 0, H],|l6j]), we use the following zenith angle (9 Z ) distributions of leptons: sub-GeV 
e-like and //-like events (SGe and SG//), divided in 10+10 bins; multi-GeV e-like and //- 
like events (MGe and MG//), divided in 10+10 bins; upward stopping and through-going 
// events (US// and UT/i), divided in 5+10 bins. The calculation of the theoretical events 
rates i?^ heo in each of the 55 bins is done as in |2l(. The SK statistical analysis is 

considerably improved with respect to 

0,1 

[21] , where only three sources of systematic errors 
were used (leading to correlated uncertainties in energy, angle, and flavor mis-identification 
+9j). Now the set of systematic errors has been enlarged to 11 entries, leading to a more 
complex structure of correlated errors affecting the i?^ heo, s, with allowance for various kinds 
of (relative) normalization and shape uncertainties. As emphasized in [22j for the solar 
neutrino case, systematic uncertainties can be implemented in the x 2 statistics through 
two equivalent methods: 1) by building the total covariance matrix ("covariance method"), 
or 2) by adding quadratic penalties in the systematic pulls ("pull method"). The latter 
approach (adopted in this work) allows to study how systematic errors alter the theoretical 

predictions from the central values i?J; heo to "shifted" values i?|f e °, in order to match the 
data. The difference R n — i?^ heo is thus useful to gauge the size and the direction of 
systematic effects in the data fit. More precise definitions and technical details are given in 
Appendix A. 

Concerning the K2K data, we use the absolute spectrum of muon events in terms of the 
reconstructed neutrino energy E 0,0, which provides a total of 29 events (here divided in 
6 bins). In this sample, the parent neutrino interactions are dominantly quasi-elastic (QE), 
and the reconstructed energy E is thus closely correlated with the true neutrino energy E u . 2 
In each bin, we attach to the theoretical prediction A^; heo the two leading systematic errors 
(due to near-far and normalization uncertainties) through the pull approach. As for the 

SK analysis, we also discuss the systematically shifted theoretical predictions N n . The 
relevant technical details are given in Appendix B. 



III. STANDARD OSCILLATIONS 

In this section we discuss updated bounds on the parameters (Am 2 , sin 2 29), governing 
the scenario of standard oscillations. The most important result is a significant strengthening 
of the upper bound on Am 2 induced by K2K data. 

Figure shows the joint bounds on the (Am 2 , sin 2 29) parameters from our analysis 
of SK, K2K, and SK+K2K data, derived through A% 2 cuts around the \ 2 minima. The 
minima and their positions are reported in the "standard oscillation" columns of Table 
The global SK best fit (x 2 min = 34.8) might seem somewhat too good, as compared with the 
A+jf = 55 — 2 degrees of freedom. Howewer, this value is only 1.8a below the typical \ 2 
expectations (A+jf ± \/2 A DF ), and is thus not suspicious from a statistical viewpoint. 

In Fig. JTJ the bounds in the left panel are very close to the official SK ones, as presented 
in 0, 13. Il6|. The bounds in the middle panel are instead slightly weaker than the official 
K2K ones [8J], especially in terms of sin 2 29. In particular, we do not find a lower bound 



2 We cannot use the full K2K data sample (56 events |8|]), which contains 27 additional events of dominantly 
non-QE origin, whose analysis is basically not reproducible outside the K2K collaboration. 
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TABLE I: Coordinates and values of the absolute \ 2 minima for separate and combined fits to SK and 
K2K data, in the scenarios of standard oscillation, pure decoherence, and oscillations plus decoherence. 



Best fits 


Standard oscillations 
SK K2K SK+K2K 


Pure decoherence 
SK K2K SK+K2K 


Oscillations plus decoherence 
SK K2K SK+K2K 


Am 2 (eV 2 ) 
sin 2 29 
p? (cV 2 ) 

x 2 


2.72xl0" 3 2.65X10" 3 2.63xl0" 3 
1.00 0.82 1.00 

34.8 9.4 45.2 


1.00 0.86 1.00 
2.42X10" 3 2.46X10" 3 2.44xl0~ 3 
42.1 11.5 54.1 


2.72xl0" 3 2.65X10" 3 2.63xl0" 3 
1.00 0.82 1.00 
0.00 0.00 0.00 
34.8 9.4 45.2 



on sin 2 29 at 99% C.L. (for ]Vdf = 2). The reason is that we cannot use the additional 
(dominantly) non-QE event sample of K2K (27 events), which would help to constrain the 
overall rate normalization and thus sin 2 29. This fact might also explain why we find the 
K2K best fit at sin 2 29 = 0.82 rather than at 1.00 as in |8|. 3 Due to the slight anticorrelation 
between Am 2 and sin 2 29 in the K2K analysis (middle panel of Fig. Q), an increase of sin 2 29 
from 0.82 to 1.00 generates a slight decrease in the K2K favored range for Am 2 , which is 
reflected in the SK+K2K best-fit Am 2 (see Table HJ). 

By comparing the left and right panels of Fig. ^ the main effect of K2K appears to be the 
strengthening of the upper bound on Am 2 , consistently with the trend of the first K2K data 
(rate only 0, no spectrum) [ijj. The main reason is that, for Am 2 ~ (4-6) x 10~ 3 eV 2 , the 
first oscillation minimum would be located at — or just above — the K2K energy spectrum 
peak, implying a strong local and overall suppression of the expected events, contrary to the 
K2K observations. 

Figure 121 shows the SK and SK+K2K bounds on Am 2 , when the sin 2 29 parameter is 
projected (minimized) away. The linear scale in Am 2 makes the K2K impact on the up- 
per limit more evident. Notice that, up to ~ 3cr, the global (SK+K2K) \ 2 function is 
approximately parabolic in the linear variable Am 2 , so that one can properly defined a 
one-standard- deviation error for this parameter. This useful feature of the SK+K2K fit 
was already argued on the basis of a graphical reduction of the official SK and K2K likeli- 
hood functions [23|, and is here confirmed through a full analysis. By keeping only the first 
significant figure in the error estimate, a parabolic fit provides the ±lcr range 

Am 2 = (2.6 ± 0.4) x 10" 3 eV 2 , (6) 

with ±Na errors scaling linearly with N (up to N ~ 3). 4 

The bounds on sin 2 29 are instead entirely dominated by SK. Figure H3 shows the A% 2 
function in terms of sin 2 29, for Am 2 projected (minimized) away in the SK fit. In this figure, 
the addition of K2K data would insignificantly change the bounds (not shown), which thus 
hold for both the SK and the SK+K2K fit. Also in this case, the nearly parabolic behavior 
of Ax 2 allows to properly define a la range, 

sin 2 2^ = 1.00^, (7) 



3 In any case, this difference is not statistically significant, since the x 2 increase from our K2K best-fit point 
to maximal mixing is < 1. 

4 This range is consistent with the estimate Am 2 — (2.7 ± 0.4) x 10 -3 eV 2 of [23j, derived from a less 
detailed analysis. 
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with the lower Na error scaling linearly with N (up to iV ~ 3). 5 Equations © and (J7J) 
concisely review the current fit to the standard oscillation parameters, as anticipated in the 
Introduction. 

Figure 0] shows the comparison between observations and best-fit predictions for the SK 
zenith distributions. The best fit point refers to the case of SK+K2K standard oscillations 
in Table [I] Since the very good agreement between data and theory is no longer a surprise, 
in the following we comment on the "fine structure" of the SK data fit. This requires, 
however, that the reader can grasp the difference between theoretical predictions with and 
without the shifts induced by correlated systematics (solid and dashed histograms in Fig. 
respectively), which is explained in detail in Appendix A. 

In Fig. 0J the comparison between solid and dashed histograms shows that the systematic 
shifts are often comparable in size to the statistical errors, implying that just increasing the 
SK atmospheric v statistics will hardly bring decisive new information on the standard 
oscillation scenario (or on physics beyond it). In the SG and MG samples, the fit clearly 
exploits the systematic uncertainties to increase the e-like event normalization, especially in 
the upward direction, so as to reduce the "electron excess" possibly indicated by SK data. 6 
Concerning /i-like events in the SG and MG samples, the fit show an opposite tendency to 
slightly decrease the normalization of (especially down-going) events. The tendency appears 
to be reversed in the high-energy UT sample. Taken together, these opposite shifts of e- 
like and /i-like expectations in the SG and MG samples seem to suggest some systematic 
deviation from the predicted fi/e flavor ratio which, although not statistically alarming (see 
the last paragraph of Appendix A), should be kept in mind. In fact, deviations of similar 
size might have their origin in neutrino physics beyond 2z/ oscillations, e.g.: 1) subleading 3z/ 
oscillations, which could slightly increase the e-like event rates (see, e.g., [T^J 2 ^ 291) . or 2) 
subleading (short-baseline) — > oscillations in the so-called (and still alive UEl) 3+1 



scenario, which mi ght slightly alter the down-going muon rates, as well as the front-detector 
K2K rate [32[ l33l. l34|. Unfortunately, since such effects are typically not larger than the 
systematic shifts in Fig. 0J they are likely (if any) to remain hidden in higher-statistics SK 
data, unless a significant reduction of the current systematics can be accomplished. The 
happy side of the story is that, for the same reasons, typical subleading effects beyond 
standard 2v oscillations do not significantly alter the fit results in Eqs. (JBJ) and (J7J). 

Figure shows the comparison between data and theory for the K2K absolute spectrum 
of events, for the same oscillation best-fit point as in Fig. 0] In this case, the amount of 
systematic deviations preferred by the fit is much smaller than the current statistical error, 
implying that there is a great potential for improvements with higher K2K statistics. 



IV. PURE DECOHERENCE (AND COMPARISON WITH OSCILLATIONS) 

In this section we discuss the bounds on the parameters (/i 2 , sin 2 29), governing the 
scenario of pure decoherence [Eq. Ej. We also discuss how this scenario compares with 
standard oscillations in fitting SK and K2K data. 



5 The "upper" (null) error is, of course, trivial. 

6 This tendency could be exacerbated by adopting the most recent at mosp heric flux calculations, which con- 
sistently predict a lower normalization at relatively low energies j24l l2fJ. l2fil l27T | . Nevertheless, preliminary 
results indicate that the induced changes in the (Am 2 , sin 2 28) bounds would be minimal 28]. 
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Figure IH1 shows the joint bounds on the (/x 2 , sin 2 29) parameters from our analysis of SK, 
K2K, and SK+K2K data, derived through Ax 2 cuts around the x 2 minima, as reported 
in the "pure decoherence" columns of Table |U Such bounds are clearly dominated by SK. 
The best fit values of /z 2 are numerically close to the previous ones for Am 2 , since for both 
oscillations and decoherence one gets large disappearance effects for xhjE v ~ 0(1), where 
x = fi 2 or Am 2 jl~D |. 

The global SK+K2K best fit for pure decoherence (Xmin = 54.1) is ~9 units higher than 
for standard oscillations (Xmm = 45.2), where ~7 units are provided by SK alone. 7 On the 
one hand, this difference shows that the data start to have some sensitivity to decoherence 
effects, and tend to globally disfavor them. On the other hand, this sensitivity is not sufficient 
to claim rejection of such effects on a purely phenomenological basis. The features that tend 
to disfavor decoherence are better explained in terms of the best-fit theoretical distributions. 

Figure [7| shows the SK zenith distributions for the global (SK+K2K) best fit to pure 
decoherence (see Table HJ), and should be compared with the analogous Fig. E|for the case 
of standard oscillations. The differences are hard to detect at first sight. In the decoherence 
case, the only noticeable difference is a slightly less pronounced suppression of the muon rates 
for increasing L, which leads to slightly less tilted \i distributions, as compared with standard 
oscillations. Notice also that the pure decoherence formula (0) implies P^ c > 1/2 always, 
while the standard oscillation formula (J1J admits minima as low as P°^ c = 0, which can 
help to get a more efficient suppression (unless the limit of averaged oscillations is reached). 
The small differences in the upgoing muon distribution slopes are then responsible for the 
Ax 2 — 7 difference between the two scenarios in SK. We do not expect that additional NC- 
enriched data in SK can add significant contributions to this difference, since (contrary to 
the decay scenario [l7j ) decoherence preserves the number of active neutrinos, as previously 
observed. Moreover, given the size of the systematic shifts in both Figs. [7| and 01 it seems to 
us that higher SK statistics will not be decisive to disentangle decoherence from oscillations. 
We conclude that SK tend to disfavor the pure decoherence hypothesis through its accurate 
measurements of the upgoing muon distribution shapes; however, there seems to be little 
room for significant improvements with longer SK exposure. 

Figure |H1 shows the K2K absolute spectrum of events for the global (SK+K2K) best fit 
to pure decoherence (see Table EJ), to be compared with the analogous Fig. EJfor standard 
oscillations. The differences between the theoretical predictions in these two Figures are 
significant, although the large statistical error bars reduce their effect to a mere A% 2 ~ 2 
difference between the oscillation and decoherence fits. The inequality P^ c > 1/2 is crucial 
to provide a slightly worse fit in the decoherence case, since it forbids a significant overall 
suppression, which can instead be more easily achieved with oscillations, particularly at low 
energy. Therefore, there is room for significant improvements in the discrimination of the 
two scenarios with higher K2K statistics, by looking at the low-energy part of the spectrum; 
in fact, a local rate suppression < 1/2 would definitely rule out the decoherence hypothesis. 

Higher K2K statistics will also allow a less coarse binning in the spectrum analysis. For 
this reason, we show in Figs. El and El the unbinned K2K spectrum for representative cases 
of standard oscillations and pure decoherence, respectively (at maximal mixing). The two 
scenarios show increasing differences for decreasing energy in all cases. Moreover, within 
the standard scenario (Fog. EJ), the spectrum shows rapid shape variation with Am 2 . A 



7 The SK collaboration finds a larger difference (A% 2 = 10.5 \4. Il5| ') by using more data. Notice that 
pure decoherence satisfies the goodness-of-fit test [Xmin/^DF ~ 0(1)] in both fits (SK and ours). 
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judicious choice of binning, especially in the low-energy part of the spectrum, might thus 
enhance the K2K discrimination power in future data analyses, and might also improve the 
determination of Am 2 in the standard oscillation case. 



V. OSCILLATION PLUS DECOHERENCE 

In this section we consider the general case of oscillations plus decoherence [Eq. (J3J)], in 
order to check whether subdominant (rather than dominant) decoherence effects can help 
the fit to the data. It turns out that, by leaving (Am 2 , sin 2 29, /i 2 ) free, the best fit is 
reached for /x 2 = for both SK, K2K, and SK+K2K, as reported in the last three columns 
of Table HJ The corresponding joint bounds on the model parameters are shown in Fig. [TT1 
through projections onto the three coordinate planes. Standard oscillations (fi 2 = 0) and 
nearly maximal mixing are clearly favored. However, the pure decoherence limit (Am 2 = 0) 
is still marginally allowed by the data. As discussed in the previous section, further K2K 
statistics might help to disfavor this limit with higher confidence. 

Since /i 2 = is consistently preferred, it makes sense to derive upper bounds on /x 2 , 
by projecting away the other two parameters (Am 2 , sin 2 29) in the fit. Figure H21 shows 
the results of this excercise, in terms of the function Ax 2 (/i 2 ). The upper bound at 3cr 
corresponds to fi 2 ~ 3 x 1CT 3 eV 2 . Values of /i 2 ~ few x 1CT 3 eV 2 are sufficiently low to 
be compatible (as remarked in with the non-observation of — > v T appearance in 
CHORUS and NOMAD @. 

We conclude by observing that a really general "oscillation plus decoherence" scenario 
should be performed with (at least) three neutrino families, in order to incorporate also solar 
and reactor neutrino flavor transitions, and possibly other phenomenological constraints. 
Models of this kind would imply several new decoherence parameters, as well as nontrivial 
complications related to unavoidable matter effects. However, given the robust theoretical 
and phenomenological basis for standard oscillations, and the lack of signals in favor of 
decoherence effects (so far), a time-consuming analysis of such extended models is perhaps 
unwarranted at present. It could become interesting, however, if further K2K data, plus 
future long-baseline accelerator data, would fail to show a clear oscillation pattern. 



VI. SUMMARY AND CONCLUSIONS 

We have analyzed in detail the current SK atmospheric neutrino data and and the first 
K2K spectral data, in order to review the status of standard — > v T oscillations, as well 
as of a possible "rival" scenario based on neutrino decoherence. We have provided updated 
bounds for the standard oscillation parameters [Eqs. (JTJ) and (J2J)], and have found no evidence 
for decoherence effects. However, the SK+K2K data are not accurate enough to rule out the 
pure decoherence case yet. The (currently weak) indications against such effects appear to 
be currently dominated by SK data, but the statistical analysis of the uncertainties reveals 
that K2K will lead further progress in this field, especially through higher-statistics tests 
of the low-energy spectrum bins. In conclusion, standard — > v T oscillations are in good 
health, and seemingly no alternative scenario is able to provide a better fit to the SK+K2K 
data. However, the "second best" explanation of SK+K2K data (i.e., pure decoherence) is 
still statistically acceptable, and can provide a useful phenomenological benchmark to test 
the emergence of the — still hidden — standard oscillation pattern. 
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APPENDIX A: STATISTICAL ANALYSIS OF SK DATA 



In this paper, the calculation of the SK zenith distributions of e-like and //-like event 
rates is performed as described in previous works 2, 2^, 21 . In the absence of oscillations, 
we normalize our absolute rates to the official 92 kTy SK estimates 0], which are currently 
based on the one- dimensional v fluxes by Honda et al. [36| for all classes of events. 8 

Our SK statistical analysis, previously based on a simplified approach to error correlations 
fl9| . has been significantly improved in two aspects: 1) it has been made more consistent 
with the error estimates of the SK collaboration, as reported in j3^| (see also jl^, Eel lilf 
and 2) it has been cast in a "pull" rather than "covariance" form, as advocated in |22| for 
solar neutrinos. 

We remind that the "pull" approach to correlated systematic uncertainties (see 2^1 and 
references therein) amounts to shift the theoretical rate i?J; heo in the n-th bin through a set 
of deviations £fcc£, where c£ is the la error associated to the k-th source of systematics (in 
the present case, k = 1, . . . , 11), and the £&'s are random variables, 



11 



ntheo 



k=l 



(Al) 



The x 2 function is then obtained by minimization over the ^'s, 



Xsk — mm 



55 
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pthco Dcxpt \ 

E ( n * n 



n=l 



a 



k,h=l 



(A2) 



where cr^ tat represent the statistical errors (or, in general, the quadratic sum of all uncorre- 
lated errors) for the n-th bin. Notice that, with respect to the discussion in j^, we now 
give allowance for correlations phk of the systematic error sources. 9 The above \ 2 function 
is exactly equivalent to build a quadratic form in the unshifted differences (-Rjf 60 — -R^ xpt ), 
with a covariance matrix defined as 



£ stat stat _j_ 



11 

E PhkC h n c k m 

k,h=l 



(A3) 



The minimization in Eq. (|A2|) leads to a solvable set of linear equations in the ^'s, whose 
solution £ fc can provide useful information on the role of systematic uncertainties in the fit 



The one-dimensional v fluxes by the Bartol group are used by the SK collaboration for comparison 
withjjhe default Honda fluxes. 
9 In 22j all error sources were strictly independent, so that phk — 5hk- Generalizing to phk 7^ ^/ifc leads to 
the appearance of the inverse correlation matrix [p -1 ]^ in the last term of Eq. 1|A2(1 . 
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22l | . For instance, it can be interesting to study how much the best-fit shifted rates 



?c=i 

differ from the standard ones (i?^ heo ), as a consequence of the systematic pulls £ fc . 10 This 
valuable information would be lost in the covariance approach. 

Our characterization of the 11 sources of SK uncertainties follows (or, at least, is inspired 
by) the detailed error description given in Q (see also 3^, 4(1 4l[ ) - There is, however, one 



important difference, related to the fact that the SK analysis is based on Monte Carlo (MC) 
simulations, while ours is based on direct calculations of the lepton rates. More precisely, 
while in the SK official anal ysis the systematic uncertainties act as re-weighting factors for 
each MC event l4l| . in our case they act directly on the observable rates R n 

(n = 1, . . • , 55), and are thus completely characterized by the 55 x 11 matrix 



Table ITT1 shows our numerical assignments for the c^'s (basically derived from |28j), inter- 
preted as percentage errors of the f?J; heo values (with or without oscillations). The second row 
of Table ITTlreports the notation of |38J for the same error sources, which we now discuss. The 
1st systematic error source (a) represents a 25% overall normalization uncertainty, due to at 
least two components: atmospheric v flux normalization error (~20%) |25[ and v cross sec- 
tion uncertainties J~ 15%), added in quadrature j^. The 2nd error source (5) is the "slope" 
uncertainty (5% [38() in the atmospheric v energy spectrum, leading to a 10 05 ~ 1 + 12% 
normalization change for a lcr "tilt" of the spectrum slope over one energy decade. Since the 
SG and UT parent neutrino spectra are separated by about two decades in energy, we take 
for them (? n = —12% and +12%, respectively. The MG and US event samples have roughly 
intermediate energies between the SG and UT ones (on a log-E^ scale), so we assign them 
c 2 n = —4% and +4%, respectively. The 3rd and 4th error sources (/3 S and (3 m ) are the overall 
fi/e flavor ratio uncertainties for the SG and MG samples, respectively. The relatively large 
ring-counting error in the MG sample makes (3 m greater than /3 S j38|. The 5th error source 
(p) is the relative normalization error between partially and fully contained (PC and FC) 
events, estimated as ~10.5% [38j]. Since only the MG/i sample contains PC events (~ 57%), 
the total (FC+PC) change in the MG/i normalization is reduced to 0.57 x 10.5% ~ 6%, as 
indicated in Table ITT1 The 6th error source (p s ) is the normalization uncertainty of upgoing 
muon events (US and UT) with respect to lower-energy events, mainly due to the relative 
cross section uncertainties. The 7th error source (pt) gives further allowance for a normal- 
ization shift of the US event rate only, dominated by the track length cut uncertainty. The 
8th and 9th error sources (rj s and i] m ) allow a small up-down asymmetry (about ±3% at la) 
in the SG and MG zenith distributions, due to geomagnetic and instrumental uncertainties. 
The 10th error source (k/) allows a small horizontal/vertical ratio uncertainty (about ±2%) 
in the zenith distributions of FC events with momentum > 0.4 GeV for both the SG and 
MG samples, mainly related to atmospheric v flux calculation uncertainties [38]. This error 
is reduced by a factor of ~ 2 in the SGe and SG/i samples, due to a ~ 1/2 component of 
low-momentum FC events which are unaffected by this uncertainty, since they loose memory 
of the neutrino direction. A similar Kf error reduction applies to the MG/i sample, due to 



10 We observe that the SK collaboration basically shows the shifted i?!^°'s in the best-fit plots, while we 



have graphically reported the unshifted i?^ heo 's in previous works |19ll20ll21| . For the sake of comparison, 
Fig. 0| and \7\ in this paper show both representations. 



10 



TABLE II: Set of lc systematic errors {c n } generated by several sources (k = 1, . . . , 11) and affecting in different ways the 
binned rates of the SK zenith distribution (n = 1, . . . , 55). All errors are given as percentage fraction of the theoretical rate in 
each bin: s„ = 100 X c*/i?^ heo . The formal identification with the uncertainties (a, 5, . . . , ««) of l38l is given in the 2nd row. 

Sec the text for further details. 
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the substantial presence of (unaffected) PC events The 11th error source (k u ) gives 

allowance for a horizontal/vertical ratio uncertainty in the UT/i sample, which is instead 
argued to be irrelevant in the US/i sample (see jl^)- Finally, we remark that the SK analysis 



contains a further error source (so-called L/E v uncertainty, ~15% |4fJ, |4l| ) that we 

do not implement in our analysis. In particular, we think that the spread in L should be 
more properly included through integration over the known neutrino pathlength distribution 
(i.e., over the neutrino production heights at different angles) rather than being interpreted 
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TABLE III: Correlation matrix phk of systematic error sources (h,k = 1, . . . , 11) in the SK analysis, as 
taken from j98| . 
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TABLE IV: Values of the SK systematic pulls {£/,} at the 


global (SK 


+K2K) best-fit point for standard 




— > v T oscillations. 
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12 3 4 
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11 


Ik 


+0.19 +0.22 -1.05 -0.99 


-0.53 


-0.03 


-0.12 


-0.99 


-0.82 


+0.07 


+0.01 



as an additional uncertainty. 11 We also observe that, in general, a U L/E V error" becomes 
meaningless in scenarios where P w is not a function of L/E u , but of L and E v separately 
(as is the case, e.g., in the presence of matter effects). 

The 11 systematic error sources listed in Table ITT1 are not necessarily independent. For 
instance, the horizontal/vertical ratio uncertainties i] s and t] m are largely due to common 
geomagnetic and instrumental effects, and are thus expected to be positively correlated. On 
the other hand, cross section uncertainties (affecting in different ways quasi-elastic and deep- 
inelastic scattering events) can generate a "migration" of events from one class to another, 
thus inducing negative correlations among the corresponding normalization uncertainties. 

Table InTl reports explicitly the correlation matrix phk among the SK systematics, as taken 
from the careful estimates in [38|. This matrix enters in the x|k definition of Eq. (jA2j) . 

Finally, Table HVI reports the systematic SK pulls {^}fc=i,...,n at the global (SK+K2K) 
best fit for standard z/ M — > v T oscillations, which have been used to calculate the "shifted" 
theoretical rates in Fig. |3J They appear to be relatively small for any k. Therefore, the 
induced differences R n — i?J; heo between the solid and dashed histograms in Fig. |U represent 
tolerable systematic shifts of "typical" size. 

APPENDIX B: CALCULATION AND STATISTICAL ANALYSIS OF THE K2K 
SPECTRUM 

In this appendix we describe our calculation of the K2K spectrum of events and our 
approach to the K2K statistical analysis. 



We understand that future SK atmospheric v data analyses might indeed fully incorporate the spread of 
L, as derived from recent three-dimensional v flux calculations (T. Kajita, private communication). 
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TABLE V: Set of systematic errors {c^} generated by several sources (k — 1, ...,7), and affecting in 
different ways the binned rates of the K2K energy distribution (n — 1,...,6). All errors are given as 
percentage fraction of the theoretical number of events N^ 100 in each bin: = 100 x c^/N^ c ° . The energy 
interval and the experimental number of events N° xpt are given in the second and third row, respectively. 
See the text for details. 
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1. Calculation of the K2K spectrum 

Our oscillation analysis of the K2K spectral data is based on 29 single-ring /i-like (lR/i) 
events 0, binned into 6 intervals of the reconstructed neutrino energy E. The experimental 
number of events A^ xpt in each bin is given in Table 13 together with the corresponding 
systematic uncertainties (discussed below). The no-oscillation expectations correspond to a 
total of 42.5 events, as graphically derived from Fig. 2 in |8j,ll6[. This data sample contains 
mainly muons from quasi-elastic (QE) scattering + n — > \i + p. 12 

In principle, for QE scattering on neutrons at rest, the reconstructed neutrino energy E 
depends only on the muon kinematics 



m n E,. —mf,/ 2 

E = ^ ^ . (Bl) 

m n - E^ + p^ cos 6» M 

However, a number of effects can produce significant deviations of E from the true neutrino 
ener gy E„, including nuclear effects (Fermi motion, Pauli blocking, proton re-scattering 
[lii I43L 0|), detection uncertainties (SK resolution in energy and angle [lH and 
contamination from non-QE events [47, 48] . Therefore, the reconstructed v energy E is 
distributed around the true v energy E v through a certain (normalized) probability distri- 
bution function R(E,E U ). This function smears out the oscillation pattern and must thus 
be taken into account, at least approximately, in the K2K analysis. 

Although the function R(E, E v ) is not explicitly provided by the K2K Collaboration 
its main fixtures can be roughly recovered from other studies. In particular. i^E, E,,) has 



been derived (for a SK-like detector) through numerical simulations in |42j, |45(. It turns 
out that jiH Hij]: 1) E is typically biased towards slightly lower values than E v (with 
(E„ — E) = b ~ few x 10 MeV), 13 and 2) the r.m.s. difference (E — E^)^ = <je increases 
from zero (at E ~ 0.2 MeV) up to a plateau (Je/E ~ 20% 0, 13 . We approximately 
embed these features in a Gaussian p.d.f. of the kind 



1 _\( E-E v +b \' 

R(E,E V ) ~ e 2 \ °* J 



(B2) 



12 We do not use the additional 27 events forming the total K2K data sample (56 events), since the simulation 
of (dominantly) non-QE events is not possible without very detailed K2K experi men tal information. 

13 Contamination from non-QE events also produces a bias b in the same direction |47| . 



13 



where we take b ~ 0.05 GeV for the bias parameter and 

'0.2 -E u [GeV]' 



a E /E v ~ 0.2 



1 — exp 



0.08 



(B3) 



for the energy reconstruction error a E above E v > 0.2 GeV (an = otherwise). Although 
we do not include the more detailed features investigated in |42j, |45J] (e.g., asymmetric tails), 
the above function R(E, E v ) appears to be good enough for our current analysis. A better 
treatment will be possible when a description of R(E, E v ) (or of equivalent information) in 
K2K will be made publicly available. 14 

The .E'-spectrum of events expected in the presence of oscillations is calculated as 



dN' 



thco 



dE 



J dE v S(E v )R(E,E u )P^(E v 



(B4) 



where Si^Ey) is the unoscillated v spectrum at the SK detector (in terms of the true neutrino 
energy). We find that, in the relevant energy range E v ~ [0.2, 3.25] GeV, the K2K spectrum 
S(E U ) is well approximated by the simple function 



S{E V ) ~ 7 (£, - 0.2) Q (3.25 - 



(B5) 



where a = 1.0056, (3 = 3.144, and E v is given in GeV. The factor 7 = 0.7389 is then 
fixed by the normalization constraint / dE v dE S(E U ) R(E, E v ) = 42.5 for no oscillation. 
The reader is referred to Fig. |U] for representative theoretical spectra dN theo /dE, with and 
without oscillations (at maximal mixing). Such spectra are in reasonable agreement with 
the corresponding K2K ones as shown in HT 



2. Statistical analysis of the K2K spectrum 

We bin both the K2K theoretical spectrum dN theo /dE and the experimental spectrum 
into six energy intervals (see Table IVj). The statistical analysis of the K2K absolute number of 
events {-/V^ } n=1 _ 6 (X = theo, expt) is performed, in analogy with the previous Appendix, 
by shifting the predictions through a set of (7, in this case) systematic deviations 



theo 



N, 



theo 



fc=i 



(B6) 



which are then minimized away in the \ 2 . The only formal difference with respect to the 
SK data fit is the use of Poisson statistics, required to deal with small number of events: 



Xk2K = mi n 



2EK 



theo 



jyexpt 



Artheo 

iV exp In n 



n=l 



n: 



expt 




\hk 



(B7) 



14 We also note that several issues related to the characterization of the E — E v difference in QE and non-QE 
events are currently being re-examined, in the light of the accuracy needed to analyze hi ghe r-statistics 
data in the K2K and future low-energy accelerator neutrino experiments (see, e.g., |44 l48l I49j1. 
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TABLE VI: Correlation matrix p^k {h, k = 1, . . . , 7) of systematics in the K2K analysis. The nontrivial 
6x6 block, related to near- far extrapolation uncertainties, is derived from [£Hj]. The 7th systematic error 
(independent overall normalization) is assumed to be uncorrelated. 
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TABLE VII: Values of the K2K systematic pulls {£ k } at the global (SK+K2K) best-fit point for standard 
v u. ~ * v t oscillations. 



k 


1 


2 


3 


4 


5 


6 


7 


Ik 


+0.02 


-0.01 


+0.29 


-0.70 


+0.38 


+0.24 


+0.38 



The minimization in Eq. (|B7|) . which can be easily performed through linearization in the 
small parameters fixes the best-fit systematic pulls The shifted theoretical predictions 
for the number of events in each bin read then 

K e ° = N^° + j:^c k n . (B8) 
fc=i 

In the systematic error budget, we include the following (dominant) error sources, for 
which one can find a detailed description in The near-far extrapolation procedure in 
K2K leads to six systematics errors (k = 1, . . . , 6) in each bin, with significant bin-to-bin 
correlations (HjJ. The overall normalization uncertainty adds a further 5% error {k = 7) 
with full correlation in each bin. The numerical values of the la systematics c£ and of the 



correlation matrix phk, as derived from |5Jj, are reported in Tables fVl and IVI[ respectively. 
These values enter in the Xk2K definition of Eq. (jB7|) . 

Finally, Table |VTl] reports the K2K systematic pulls {C, k }k=i,...j at the global (SK+K2K) 
best fit for standard — > v T oscillations, from which the "shifted" theoretical rates in Fig. |3] 
are calculated. All the pulls appear to be relatively small. 
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FIG. 1: Standard oscillations in the — > v T channel: Bounds on the parameters (Am 2 , sin 2 29) 
from SK atmospheric data (left panel), K2K spectral data (middle panel), and their combination 
(right panel). The solid and dotted curves refer, respectively, to the 90% and 99% C.L. contours 
for iV DF = 2 (Ax 2 = 4.61 and 9.21). 
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FIG. 2: Standard oscillations in the v„ — > z/,- channel: Bounds on Am 2 for unconstrained sin 2 29 
from SK (dashed curve) and SK+K2K (solid curve). The intersections with the horizontal dotted 
lines give the 2a and 3a (upper and lower) bounds on Am 2 for Ndf = 1- By fitting the SK+K2K 
curve with a parabola, the ±la interval is derived as Am 2 = (2.6 ± 0.4) x 10~ 3 eV 2 . 
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FIG. 3: Standard oscillations in the — > z/,- channel: Bounds on sin 2 26 for unconstrained Am 2 
from SK data. The inclusion of K2K data induces negligible changes (not shown). The intersections 
with the horizontal dotted lines give the 2a and 3a lower bounds on sin 2 26 for iVrjF = 1- By fitting 
the curve with a parabola, the ±1<t interval is derived as sin 26 = 1.00±g;gg, where only the first 
significant digit is kept in the lower error. 
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FIG. 4: Standard oscillations in the v„ — > v T channel: SK experimental zenith distributions 
(±l<7 s tat)) compared with the corresponding theoretical ones at the global (SK+K2K) best-fit 
point given in Tabled All distributions are normalized to the unoscillated predictions in each bin. 
For the theoretical event rates, we show both the central values -R^ heo (dashed histograms) and the 
"shifted" values i?^ h °° (solid histograms), which embed the effect of systematic pulls. The difference 
between R n and R^ eo shows how much (and in which direction) the correlated systematic errors 
tend to stretch the predictions in order to match the data. See the text and Appendix A for details. 
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K2K energy spectrum, standard oscillations 
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FIG. 5: Standard oscillations in the — > v T channel: Absolute spectrum of (dominantly QE) 
events in K2K, as a function of the reconstructed neutrino energy E. The data points (29 events 
total) are shown as dots with ±lcr s t a t in each of the six bins. The dotted histogram represents 
our calculations for no oscillation. The solid and dashed histograms represent the theoretical 
predictions iV* heo and iV* heo at the global (SK+K2K) best-fit point in Table ID with and without 
systematic shifts, respectively. See the text and Appendix B for details. 
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v ^ <r> v T pure decoherence 
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FIG. 6: As in Fig. ^ but for the pure decoherence scenario. Notice that the bounds on the 
(/j 2 , sin 2 29) parameters are dominated by SK. 
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FIG. 7: As in Fig.EJ but for the global (SK+K2K) best-fit point to the pure decoherence scenario, 
as given in Table [fl 
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K2K energy spectrum, pure decoherence 
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FIG. 8: As in Fig.El but for the global (SK+K2K) best-fit point to the pure decoherence scenario, 
as given in Table HJ 
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K2K energy spectrum, standard oscillations 
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FIG. 9: Standard oscillations: Unbinned K2K theoretical spectrum at maximal mixing and 
three representative values of Am 2 . 
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FIG. 10: Pure decoherence: Unbinned K2K theoretical spectrum at maximal mixing and for three 
representative values of ji 2 . 
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FIG. 11: General case (oscillations plus decoherence) : Volume allowed by SK+K2K in the param- 
eter space (Am 2 , sin 2 29, fj 2 ) at 90% and 99% C.L. for N DF = 3 (Ax 2 = 6.25 and 11.34), shown 
through its projections onto the coordinate planes. The best-fit is reached in the limit of standard 
oscillations (fi 2 = 0). However, the case of pure decoherence (Am 2 = 0) is still marginally allowed. 
In all cases, the mixing is nearly maximal (sin 2 26 ~ 1) 
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Bounds on \x for unconstrained (Am 2 , sin 2 2i^) 
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FIG. 12: General case (oscillations plus decoherence) : Upper bounds on the decoherence parameter 
/i 2 , for unconstrained (5m 2 , sin 2 29). The intersections with the horizontal dotted lines give the 2a 
and 3a upper bounds on [i 2 for Nbf = 1- 
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